function [KxR1 KxR2] = funcKxR(nx, nz, zeta, CONSTS)

    [KxS1 KxS2] = funcKxS(nx, nz, zeta, CONSTS);

    Kx = funcKxfull(nx, nz, zeta, CONSTS);
    
    KxR1 = Kx - KxS1;
    KxR2 = Kx - KxS2;
    
end